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ABSTRACT 

This paper is concerned with a new development of a 
delay-robust stabilizing feedback control law for 
linear ordinary differential equation coupled with two 
linear first order hyperbolic equations in the actuation 
path. A second change of variables that reduces the 
stabilization problem of the PDE-ODE system to that 
of a time-delay system for which a forecaster can be 
constructed. Hence, by choosing the pole placement 
on the ODE when constructing the forecaster, 
enabling a trade-off between convergence rate and 
delay-robustness. The proposed feedback law is 
finally proved to be robust to small delays in the 
actuation. 

Keyword: Hyperbolic partial differential equation, 
Time-delay systems, Stabilization. 

INTRODUCTION 

In this paper we develop a linear feedback control law 
that achieves delay-robust stabilization of a system of 
two hetero directional first-order hyperbolic Partial 
Differential Equations (PDEs) coupled through the 
boundary to an Ordinary Differential Equation 
(ODE). It has been observed, that for many feedback 
systems, the introduction of arbitrarily small time- 
delays in the loop may cause instability under linear 
state feedback. In particular, for coupled linear 
hyperbolic systems, recent contributions have 
highlighted the necessity of a change of paradigm in 
order to achieve delay-robust stabilization. The main 
contribution of this paper is to provide a new design 
for a state-feedback law for a PDE-ODE system that 
ensures the delay-robust stabilization. The original 
system can then be rewritten as a distributed delay 
equation for which it is possible to derive a stabilizing 
control law. 


Problem Formulation 

In this section we detail the notations used through 
this paper. For any integer p > 0, ||. || e p is the 
classical euclidean norm on R p . We denote by 
^([O.l^R), or Z. 1 ([0,1]) if no confusion arises, the 
space of real-valued functions defined on [0,1] whose 
absolute value is integrable. This space is equipped 
with the standard L 1 norm, that is, for any / e 
^([0,1]) 



We denote L 2 ([0,1],E) the space of real-valued 
square-integrable functions defined on [0,1] with the 
standard L 2 norm, i.e., for any / 6 L 2 ([0,1]), R) 

II / 11*2 = [ f 2 {x)dx. 

Jo 

The set L°°([0,1], M.) denotes the space of bounded 
real-valued functions defined on [0,l]with the 
standard U° norm, i.e., for any / 6 L°°([0,1], R) 

II / || L *= sup xe[(U] |/(x)|. 

In the following, for (u,v,X) E (L 2 ([0,1])) 2 AR p , we 
define the norm || (u,v,X) || = || u || L 2 + || v || L 2 + 
II X || RP . 

The set c p ([0,l]) stands for the space of real valued 
functions defined on [0,1] that are P times 
differentiable and whose Pth derivative is continuous. 
The set r is defined as 

T = {(*, O E [0,1] 2 S. t. f < x}. 
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c(t) stands for the space of real-valued continuous 
functions on r. For a positive real k and two reals 
a < b, a function /defined on [a, b] is said to k- 
Lipschitz if for all (x,y) G [a, b] 2 , it satisfies 
| fix') — /(y) | < k\x — y\. The symbol I P (or I if no 
confusion arises) represents the pXp identity matrix. 
We use the notation /(s) for the Laplace transform of 
a function /(t), provided it is well defined. The set <A 
stands for the convolution Banach algebra of BIBO - 
stable generalized functions in the sense of 
Vidyasagar. A function g(.) belongs to <A if it can be 
expressed as 

00 

g(x) = a A 0 + ^ 9i s(t - tj), 

i= 0 

Where g r E L 1 (E + ,E),£j> 0 \9i\ < °°, 0 = t 0 < t ± < 
••• and 5(.) is the direct distribution? The associated 
norm is 

iisiU=iiar»L*+2l*l- 

i> 0 

The set A of Laplace transforms of elements in Jl is 
also Banach algebra with associated norm 
II 9 IU=» 9 IU- 

System Under Consideration 

We consider a class of systems consisting of an ODE 
coupled to two hetero directional first-order linear 
hyperbolic systems in the actuation path. We consider 
systems of the form: 

u t (t,x) + Au x (t,x) = a + ~(x)v(t,x) 
v t (t,x) — pv x (t,x) = a~ + (x)u(t,x ) 

X(t) =AX(t) + Bv[t, 0 ), 

Evolving in {(t,x)s. 11 > 0,x G [0,1]}, with the 
boundary conditions 

u(t, 0) = qv(t, 0) + CX(t ) 
v(t, 1) = pu(t, 1) + U (t), 

Where X EMP is the ODE state, u(t,x) G 
E and v(t, x) G E are the PDE states and U (t) is the 
control input. The in-domain coupling terms a~ + and 
a + ~ belong to C° ([0,1]), the boundary coupling 
terms q =£ 0 (distal reflexion) and p (proximal 
reflexion), and the velocities A and p are constants. 
Furthermore, the velocities verify 
—p < 0 < A. 

The initial conditions of the state ( u , v) are denoted 
u Q and v 0 and are assumed to belong to Z, 1 ([0,1], E) 


and we consider only weak L 2 solutions to the system. 
The initial condition of the ODE is denoted X 0 . 
Remark that this system naturally features several 
couplings that can be extended to the case q = 0 with 
a slight modification of the backstepping 
transformation. 

Control Problem 

The goal of this paper is to design a feedback control 
law U = JC[(u,v,X)] where X: (L 2 ([0,1])) 2 XE p -> 
E is a linear operator, such that: 

The state (u,v,X) of the resulting feedback system 
exponentially converges to its zero equilibrium 
(stabilization problem), i.e. there exist k Q > 0 and 
v > 0 such that for any initial condition (u 0 , v 0 , X 0 ) E 
(L 1 [0,l]) 2 AE p 

II (u,v,X) ||< k 0 e~ vt || 

(uo,t?o»*o) II» t>0. 

The resulting feedback system is robustly stable with 
respect to small delays in the loop (delay- 
robustness), i.e. there exists 8* > 0 such that for any 
8 G [0,5*], the control law I/(t — 8) still stabilizes. 

A control law that satisfies these two constraints is 
said to delay-robustly stabilize system 

In this paper, we make the two following 
assumptions: 

Assumption 1: 

The pair (A,B) is stabilizable, i.e. there exists a 
matrix K such that A + BK is Hurwitz. 

Assumption 2: 

The proximal reflection p and the distal reflection q 
satisfy |pq| < 1. 

The first assumption is necessary for the 
stabilizability of the whole system, while the second 
assumption is required for the existence of a delay- 
robust linear feedback control. This second 
assumption is not restrictive since, if is not fulfilled, 
one could prove using arguments similar to those in 
Auriol that the open-loop transfer function has an 
infinite number of poles in the complex closed right 
half-plane. Consequently, one cannot find any linear 
state feedback law U{. ) that delay robustly stabilizes. 
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Lemma 

Take A,B and K verifying Assumption 1 and any k 
such that holds. Then, the control law U 0DE (t) = 

K B U 0DE (v)dv}, 

Exponentially stabilizes T(t). Furthermore, the state 
feedback 

UodeM = U 0DE (t ) - (p - k)CX(t - i) 
Exponentially stabilizes X(t). 

Proof 

For the state-forecaster feedback U 0DE (.), the closed- 
loop system in 

Y (t) = AY (t) + BU ode (t — , t > ^ satisfies 

F(0 = (A + BK)Y(tl t>\ 

Exponential stability is guaranteed by the fact that 
(A + BK) is Hurwitz. By construction of Y (t) and 
using 0(s) = 1 — (p — k)qe~ r , we have that X(t), 
solution of 

X(t ) - (p - /c)gA(t - r) = AX(t) - 

(p — k)qAX (t — r) + (p — 

kBCXt- t+ BUODEt- lju, satisfies for any £>r, 


X(t) = (p - k)qA(t - r) + y (t). 

Since |(p — /c)q| < 1 by \kq\ + |pp| < 1, X(t) is also 
exponentially stable. 


CONCLUSION 

In this paper, a delay-robust stabilizing feedback 
control law was developed for a coupled hyperbolic 
PDE-ODE system. The proposed method combines 
using the back stepping approach with a second 
forecaster-type feedback. The second feedback 
control is obtained after a suitable change of variables 
that reduces the stabilization problem of the PDE- 
ODE system. 
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